Quantum fingerprints of classical Ruelle-Pollicot resonances 
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Abstract 

N-disk microwave billiards, which are representative of open quantum systems, are studied 
experimentally. The transmission spectrum yields the quantum resonances which are consistent 
I with semiclassical calculations. The spectral autocorrelation of the quantum spectrum is shown to 

PsJ I be determined by the classical Ruelle-Pollicot resonances, arising from the complex eigenvalues of 

I the Perron-Frobenius operator. This work establishes a fundamental connection between quantum 

classical correlations in open systems. 

The quantum-classical correspondence for chaotic systems has been studied extensively in the 
context of universality and periodic orbit contributions. This approach has focussed on eigenvalues 
and eigenfunctions and their statistical properties. Universality has been shown to arise from Random 
Matrix Theory [||, while periodic orbit contributions have been analyzed in the semiclassical scheme 
for calculations of eigenvalue spectra |Q] and constructions of eigenfunctions . 

An entirely different approach is to consider correlations of observables. In the classical context 
a probabilistic approach is best taken with Liouvillian dynamics. In certain classical systems these 
have been shown to lead to Ruelle-Pollicot (RP) resonances 0, ^, arising from complex eigenvalues of 
the Perron-Frobenius operator. In open systems, this leads to a quantitative description of the time- 
^ . evolution of classical observables, the most common being the particle density. In the quantum context, 

\ diffusive transport has been argued to be intimately connected with Liouvillian dynamics, not just in 

' disordered systems where the correspondence is made with nonlinear cr— models of supersymmetry Q 

■ but also in individual chaotic systems which represent a ballistic limit. 

(<~^ \ In this paper we present a microwave experiment which demonstrates this deep connection between 

C3 ' quantum properties and classical diffusion. Our experiment is a microwave realization of the well- 

known n-disk geometry, which is a paradigm of an open quantum chaotic system, along with other 
systems such as the Smale horseshoe and the Baker map 0. The classical scattering function of the 
chaotic n-disk system is nondifferentiable and has a selfsimilar fractal structure. A central property is 
^ ■ the exponential decay of an initial distribution of classical particles, due to the unstable periodic orbits, 

^ . which form a cantor set, hence the name fractal repeller. The experimental transmission spectrum 

^ \ directly yields the frequencies and the widths of the low lying quantum resonances of the system 

1^, ^, which are in agreement with semiclassical periodic orbit calculations 11, The same 
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spectra are analyzed to obtain the spectral wave- vector autocorrelation C(k) |Q. The wave vector 
dependence of the spectral autocorrelation is shown to be completely described by the leading RP 
resonances of the corresponding classical system. The small k (long time) behavior of the spectral 
autocorrelation provides a measure of the quantum escape rate, and is shown to be in good agreement 
with the corresponding classical escape rate. For large k (short time), the contribution of classical RP 
resonances is observed as non-universal oscillations of the autocorrelation. Thus we are experimentally 
able to observe the classical RP resonances in a quantum experiment, for the first time. 

The experiments are carried out in thin microwave structures consisting of two highly conducting 
Cu plates spaced d ~ 6mm and about 55 x 55cm in area. Disks and bars also made of Cu and of 
thickness d are placed between the plates and in contact with them. In order to simulate an infinite 
system microwave absorber material ECCOSORB AN-77 was sandwiched between the plates at the 
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edges. Microwaves were coupled in and out using terminating coaxial lines which were inserted in 
the vicinity of the scatterers. All measurements were carried out using an HP8510B vector network 
analyzer which measured the complex transmission (S21) and reflection {Sn) S-parameters of the coax 
+ scatterer system. It is crucial to ensure that there is no spurious background scattering due to the 
finite size of the system. This was verified carefully as well as that the effects of the coupling probes 
were minimal and did not affect the results. 

In this essentially 2-D geometry, Maxwell's equation for the experimental system is identical with 
the Schrodinger time-independent wave equation (V^ + /c^)^ = with = the z-component of 
the microwave electric field. This correspondence is exact for all frequencies f < fc = c/2d = 25GHz. 
(Note that k = 27r//c, where c is the speed of light). It is this mapping which enables us to study the 
quantum properties of the 2-D systems. For all metallic objects in the 2-D space between the plates, 
Dirichlet boundary conditions apply inside the metal. Similar microwave experiments, which exploit 
this QM-E&M mapping, have been used to study quantum chaos in closed 14] and open systems 
1^]. See for details of the experiments. 

The transmission function 5*21 (/) which we measure is the response of the system to a delta- 
function excitation at point fi probed at a different point r2, and is determined by the wavefunction 
^ at the probe locations ri and f2. In our experiments the coax lines act as tunneling point contacts, 
and hence it can be shown |p that S'2i(/) = A{f)G{ri,f2,k) is just the two-point Green's function 
G{fi,r2,k). The scaling function A{f), which represents the impedance characteristics of the coax 
lines and probes, is sufficiently slowly varing and can be treated as a constant practically. Because we 
ensure that the coupling to the leads is very weak, any shifts due to the leads are negligible (< 10~^ 
of the resonance frequencies and widths)!^. The n-disk systems are investigated in the fundamental 
domain [|l^], as shown in the inset to Fig.|^, with angles 90° (n = 2), 60° (n = 3), and 45° (n = 4). 
A typical trace for the 3-disk system is shown in Fig. |l|. See for details of the comparison of the 
resonances between experiments and semiclassical calculations. 
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Figure 1: Transmission spectrum S2i{k) of the 3-disk system in the fundamental domain with disk 
separation R = 20\/3 cm and radius a = 5 cm. The dashed line is the semiclassical calculation. See 
1^] for details. (Insets) Geometry of the n-disk open billiard (n = 2,3,4). The solid lines represent 
the fundamental domain in which the present experiments were carried out. 



The spectral autocorrelation function was calculated as C{k) = (|5'2i(A; — {K/2))\'^\S2i{k + (k/2))P)^. 
The average is carried out over a band of wave vector centered at certain value ko and of width AA;. 
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Since the transmission function is the superposition of many resonances, |5'2i(A;)p = ~ 
Si)"^ + )j with Si + is[ the semiclassical resonances and Cj the couphng which depend on the location 
of the probes, we have jl^ 

C(K) = 7r> „ ,,0 , , , , „M2 - (1) 



In the case that there are no overlapping resonances, |sj — Sj| >> (s^ + s^), the small k behavior of 
the autocorrelation is C{n) ~ Tr^^j 2c?s-/(k^ + 4sf ). According to semiclassical theory, the above sum 
can be replaced by a sing le Lorentzian 0, ||] 

where 7 = 7d, the classical escape rate with the velocity scaled to 1. The above equation was used to 
fit the spectral autocorrelation for small k and thus obtain the value of the experimental escape rate 
Iqm 1^! ^ • Good agreement of the escape rate is obtained between obtained from the experiments 
and 7ci of the classical theory . 

For intermediate At, the semiclassical prediction of Eq. (|2|) fails because of the presence of the 
periodic orbits, which leads to non-universal behavior. Non-universal contributions can play in general 
a crucial role in determining the overall structure of the spectral autocorrelation, since they can be 
of the same order of the universal result of Random Matrix Theory. Recently, Agam derived a 
semiclassical theory to build the connection between the nonuniversality of the spectral autocorrelation 
and the classical RP resonances. Consider the quantum mechanical propagator [^] 

K{TuT2,t) = i^.j G{n,f2, ^2me/n^)e-'''/^de, (3) 

with e = h'^k'^ /2m. The integration is performed around £q = fi'^kQ/2m, with Ae = hvAk and 
V = hko/m is the group velocity of the classical particle. The integration in the e space can be 
changed into that in the k space as K{fi,f2,t) = (i;/27ri)e~*^o*/'* /^^, G(ri, f2, /cq + k)e~^^^''dk. The 
particle density is p{t) = \K{fi,r2,t)\^ . The autocorrelation of the particle density is Cp^r) = 
{p{t) p{t + t)) ^ — {p)1 with {p)^ = limf^oo(l/T')/o p{t)dt. Using the diagonal approximation, we 
get Cp{t) = (A/cf^/47r^y^) / dKC{n)e~™'^'^ . Here V is the volume of the system with 1/ — > 00 
for open system. If one assumes that the above correlation is classical, one has j2^, C'p(r) = 
Si^i 2&je~'''*'^'^ cos 7j'f T, where the q are the coupling coefficients, 7^ it the RP resonances of the 
corresponding classical system in wave vector space. Taking the Fourier transform of the above ex- 
pression / drC p{T)e^'^'"'^ , we get 
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C{-) = E „2 , ?ll„.2 - (4) 



with b[ = 2-nVHi/Akv^. 

We now turn to the classical dynamics of the system. The classical evolution is described by the 
Perron-Frobenius operator whose spectrum, known as the RP resonances, can be calculated as the 
poles of the classical Ruelle ^-function. For the hard disk system, the classical Ruelle C, -function is 



C/3(^) = n [1 - exp(5Lp)/|Ap|A^-i] , (5) 



here, Lp the length of the periodic orbit Ap the eigenvalue of the monodromy matrix associated 
with the periodic orbits. The C -function is analytical in the half-plane Res < — and has poles in 
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the other half-plane. In particular, Ci3{s) has a simple pole at s = —P{f3). Here, P{I3) is the so-called 
Ruelle topological pressure, from which all the characteristic quantities of classical dynamics can be 
derived in principle. The classical escape rate is 7c/ = — The poles of C^i^) with (3 = 1 are 
calculated since they contribute to the RP resonances with the sharpest width. 

For the integrable 2-disk system in the fundamental domain, there is just one prime periodic orbit. 
We have Cii^) = 1 — to, where to = exp[s(i? — 2a)]/A, and A = (cr — 1) + \Ja[o — 2), with the disk 
separation ratio a = Rja. The classical scattering resonances are s„ = (In A it i2n'K)/{R — 2a), with 
n = 1,2,- ••. The classical escape rate is 7cz = (In A)/(i? — 2a). The semiclassical resonances in 
the fundamental domain are ||2^ (2n7r + i(l/2) In A)/(i? — 2a) with n = 1, 2, • • •. Substituting the 
semiclassical resonances into Eq.(|^), one may express the full two-point correlation function as 

oo ^ 

C(k) oc -^—^ -TK, (6) 

where 7 = 7ci, 7' = 271 / [R — 2a). On the other hand, since the RP resonances of the system are 
(In A + z2n7r)/(i? — 2a), if one puts these resonances into Eq. (^, the above expression (^ follows 
immediately. The experimental RP resonances are obtained from the autocorrelation by fitting it with 
Eq. (^). Since the transmission coefficient 5*21 (/) and also the couplings Cj of the quantum resonances 
depend on the location of the two probes, so does the coupling 6„ of the classical RP resonances. The 
comparison between experiment and theory is shown as Fig. Q (left). 
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Figure 2: Autocorrelation function C(k) vs k (cm~^) and Ruelle- Pollicott resonances of the n-disk 
system in the fundamental domain, (left). 2-disk system with i? = 40 cm and a = 5 cm. (center). 
3-disk system with R/a = 4\/3 and a = 5 cm. (right). 4-disk system with R/a = 4\/2 and a = 5 cm. 
(top). The gray line is the experimental autocorrelation. The solid line is the numerical fit to Eq. 
(Q). (bottom). The dashed lines are the decomposition of the experimental C{k) into Lorentzians. 
The filled squares are the position 7^ and the bars the widths 7^- of the experimental RP resonances 
obtained from the decomposition. The open circles with dotted bars are the positions 7j and width 
7^- of the predicted RP resonances calculated from Eq. (||). The insets in the top panels show the 
leading periodic orbits (solid lines) for the different geometries. 

For the chaotic n-disk system, making use of the cycle expansion [^] and also the symmetry 
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factorization of the classical Ruelle C-function, the RP resonances can be calculated very accurately 
|24|. For the 3-disk system in the fundamental domain, the classical RP resonances are calculated by 
using 8 prime periodic orbits up to period 4. The calculation is very accurate for Res < 0.8 |2^. The 
classical RP resonances can also be obtained experimentally by fitting the autocorrelation with Eq. 
(Q. Because of the finite range of the experiment spectrum, only the first 8 or 9 RP resonances with 
small real part were obtained. The experimental autocorrelation with the theoretical prediction are 
shown as Fig. ^ (center). 

For the 4-disk system in the fundamental domain, the classical RP resonances are calculated by 
using 14 prime periodic orbits up to period 3. The experimental autocorrelation with the theoretical 
prediction are shown as Fig. ^ (right). 

The agreement between the experimental RP resonances and the theoretical ones for the 2-disk 
sytem is 6% for the positions 7j and better than 30% for the widths 7^-, is 7% for 7^ and 11% for 7^ for 
the 3-disk sytem, and is 8% for 7^ and 17% for 7^ for the 4-disk sytem. We note that these agreements, 
in particular the wave- vector locations 7^ should be considered as very good. The principal sources for 
the residual discrepancies are the nonideality of the absorbers, small symmetry-breaking perturbations 
and the suppression of some resonances at the neighborhood where the antennas are coupled which 
affects the autocorrelation function, therefore the position and the widths of RP resonances. Also very 
broad resonances are difficult to identify and can lead to an apparent enhancement of the observed 
widths, which can possibly account for the systematically larger widths that are observed. 

Our investigation clearly demonstrates that the whole spectral autocorrelation can be understood 
completely in terms of the classical RP resonances. The meaning of these RP resonances in the 
classical context can be understood as follows. If one shoots particles toward the hard disk scatterer, 
the number of particles that will remain in the scattering region will decay as N{t) = ^jOjC""'*. 
Besides the general exponential decay at uq = vjci, there are oscillations due to the fact that the RP 
resonances Oj = u(7j ± ij^) are not always real |Q as contrasted with the purely diffusive system. 
Taking the Fourier transform of N{t), one can identify the Lorentzians in the spectrum with the 
RP resonances. Our work demonstrates that suitable quantum correlations diffuse just like classical 
observables in an open system. 

It is remarkable that the same experiment yields both the quantum resonances and the classical RP 
resonances. Thus we have demonstrated experimentally the profound connection between quantum 
properties and classical diffusion. This connection is best seen in open quantum systems. While 
we have studied the model n-disk geometry, the results have broad implications for arbitrary chaotic 
geometries. The results of this work also have wider implications in a variety of phenomena in different 
fields in physics, such as photodissociation of atoms ||2^, nuclear decay [^, electronic transport, fluid 
dynamics, and acoustic and electromagnetic propagation. 

We thank P. Pradhan for useful discussions. This work was supported by NSF-PHY-9722681. 
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